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r^ Abstract 

We revisit the classical problem of nucleated polymerisation and derive a range of exact results 

^_^ describing polymerisation in systems intermediate between the well-known limiting cases of a re- 
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Oj action starting from purely soluble material and for a reaction where no new growth nuclei are 
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I. INTRODUCTION 

The classical theory of nucleated polymerisation pQ describes the growth of filamentous 
structures formed through homogeneous nucleation[2H7]. This framework was initially devel- 
oped by Oosawa and coworkers in the 1960s[ll E] to describe the formation of biofilaments, 
including actin and tubulin. This theory has been generalised to include secondary nucle- 
ation processes by Eaton and Ferrone[9] in the context of their pioneering work elucidating 
the polymerisation of sickle haemoglobin, and by WegnerjlOj in order to include fragmenta- 
tion processes into the growth model for actin filaments. 

For irreversible growth in the absence of pre-formed seed material and secondary nucle- 
ation pathways, in 1962 Oosawa presented solutions to the kinetic equations which were 
very successful in describing a variety of characteristics of the polymerisation of actin and 
tubulin. The other limiting case, namely where seed material is added at the beginning of 
the reaction and where no new growth nuclei are formed during the reaction, is also well 
known. In this paper, we present exact results which encompass all cases between these 
limiting scenarios, extending the results of Oosawa for a system dominated by primary nu- 
cleation to the case where an arbitrary concentration of pre-formed seed material is present. 
We also discuss a range of general closed form results from the Oosawa theory for the be- 
haviour of a system of biofilaments growing through primary nucleation and elongation. We 
then compare the behaviour of systems dominated by primary nucleation to results derived 
recently for systems dominated by secondary nucleation. 



II. RESULTS AND DISCUSSION 

A. Derivation of the rate laws for the polymer number and mass concentrations 

The theoretical description of the polymerisation of proteins such as actin and tubulin to 
yield functional biostructures was considered in the 1960s by Oosawa[8j. For a system that 
evolves through primary nucleation of new filaments, elongation of existing filaments, and 
depolymerisation from the filament ends, the change in concentration of filaments of size j, 



denoted f{j,t), is given by the master equation [H [8]: 
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where k+, kos, kn are rate constants describing the elongation, depolymerisation and nucle- 
ation steps and m(t) is the concentration of free monomeric protein in solution. The factor 
of 2 in Eq. ([I| originates from the assumption of growth from both ends. For the case of 
irreversible biofilament growth, the polymerisation rate dominates over the depolymerisa- 
tion rate; from Eq. ([I]), the rate of change of the number of filaments, -P(t), and the free 
monomer concentration, m{t), were shown by Oosawa under these conditions [HIH] to obey: 
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Combining Eqs. ^ and (pi) yields a differential equation for the free monomer concentration [T]: 
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Here, we integrate these equations in the general case where the initial state of the system 
can consist of any proportion of monomeric and fibrillar material; this calculation generalises 
the results presented by Oosawa to include a finite concentration of seed material present at 
the start of the reaction. Beginning with Eqs. ^ and (|3]), the substitution z{t) := log(m(t)) 
followed by multiplication through by dz/dt yields: 
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Integrating both sides results in: 
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we obtain a separable equation for dz/dt, which can be solved to yield: 
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Integration and exponentiation yields the expression for m{t): 
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Inserting the appropriate boundary conditions in terms of r?7,(0) and -P(O) fixes the values of 
the constants A and S, resulting in the final exact result for the polymer mass concentration 
M{t)=mtot-m{t): 

r / 1 M^ 

M(t) = mtot - m(0) /isech [v + XoP'^fit) (9) 

where the effective rate constant A is given by A = ^y2knk+7n{0)"■'' and /3 = 2/nc, /x = 
v/l + 72, z/ = arsinh (7) for 7 = 2A;+P(0)/(/3iA). 

We note that this expression only depends on two combinations of the microscopic rate 
constants, fco = 2A;+P(0) and A. The result reveals that A controls the aggregation resulting 
from the newly formed aggregates, whereas kg defines growth from the pre-formed seed 
structures initially present in solution. In the special case of the aggregation reaction starting 
with purely soluble proteins, -P(O) = 0, m(0) = mtot, these expressions reduce to /i — > 1 
and u -^ 0, and Eq. (|9| yields the result presented by Oosawa[l] and the single relevant 
parameter in the rate equations is A. Interestingly, generalisations of Eq. ^ which include 
secondary pathways, maintain the dependence on A and ko but introduce an additional 
parameter analogous to A for each active secondary pathway [TTHT3]. 

An expression for the evolution of the polymer number concentration, P{t) may be derived 
using Eq. ^. Direct integration of Eq. ^ gives the result for P{t): 

Pit) = P(0) + fc.m(Or /""^^'^ + ^"^y^~^""^('^^ (10) 
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Eqs. (10) and (tol) give in closed form the time evolution of the biofilament number and mass 
concentration growing through primary nucleation and filament elongation. 

B. Characteristic features of growth involving pre-formed seed material 

Insight into the early time behaviour of the polymer mass concentration can be obtained 
by expanding Eq. (|9|) for early times to yield: 



M{t) ^A M(0) + kom{0)t + m(0) [A^ - kl]ty2 + 0{t^) (11) 
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FIG. 1. Nucleated polymerisation in the presence of seed material. The thick dashed lines are 
the exact solution to the rate equations Eq. ([9]); the thin solid lines are calculated from numerical 
simulations of the master equation Eq. (IT]). The dotted lines are the initial gradients dM/dt\t=o = 
M(0) + 2A;4.?TT-(0)-P(0)t; a lag-phase exists when the initial gradient is not the maximal gradient. 



The numbers accompanying each curve are ko/X; Eq. (14) predicts that a lag-phase only exists 



when this ratio is less than unity, (a): Polymerisation in the presence of an increasing quantity 
of seed material of a fixed average length (5000 monomers per seed) added at the beginning of 
the reaction. The seed concentrations given as a fraction of the total concentration of monomer 
present are right to left): 0, 0.01, 0.04, 0.1, 0.2, 0.5. (b): Nucleated polymerisation in the presence 
of a fixed quantity (1% of total monomer in the system) of seed material of varying average length. 
The average number of monomer per seed are (right to left): N/A (unseeded), 5000, 1000, 500, 
200, 50. The other parameters for both panels are: nitot = lO/uM, ric = 3, knrn^^^ = 1 • 10~^s~^, 
k^ = l- lO^M-i s-i. 



This expression recovers the characteristic ~ t^ dependence of the Oosawa theory and has 
an additional term hnear in time relating to the growth of pre-formed aggregates. 

In many cases, Eq. ^ describes a sigmoidal function with a lag phase. The time of 
maximal growth rate, tmax, can be found from the inflection point of the sigmoid from the 
condition cPM/dt^ = 0: 
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such that a lag phase exists only for: 

artanh I W- J > arsinh (7) (13) 

Using the composition sinh(artanh(x)) = x/\/\ — x^ reduces this to the simple condition: 

A;o < A (14) 

In other words, a point of inflection exists if the growth through elongation from the ends 
of pre-existing seeds, fco? is less effective that the effective growth through nucleation and 
elongation of new material, A. This result imples that an increased nucleation rate promotes 
the existence of an inflection point, whereas an increased elongation rate or an increased 
level of seeding tends to disfavour its existence. In particular, we also note that in the 
absence of nucleation, an inflection point cannot exist in the polymer mass concentration as 



a function of time. Interestingly, the result Eq. ( 14 ) is analogous to the criterion applicable 
for fragmentation dominated growth where a lag phase only exists when the parameters 
controlling fragmentation-related secondary nucleation is larger than ek^. 
The maximal growth rate, rmax, is given by: 
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which occurs at a polymer mass concentration M^ax given from Eq.: 
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The lag time, riag := t^ax - M(tmax)/^max, is then given by: 
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Interestingly, from Eq. (17), we note that a point of inflection can never exist for P{t) 
for simple nucleated polymerisation. By contrast, when secondary pathways are active, an 
inflection point can frequently be present [12] . 
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TABLE I. Comparison of biofilament growth dominated by primary and secondary nucleation 
pathways 

C. Comparison betw^een nucleated polymerisation in the presence and absence of 
secondary pathw^ays 

Many systems that evolve through nucleated polymerisation display characteristic scaling 
behaviour [TTHT5] . This behaviour can be seen to be a consequence of the fact that under 
many conditions, the rate equations are dominated by a single parameters that corresponds 
to the dominant form of nucleation: A for classical nucleated polymerisation and k,2- for 
polymerisation in the presence of secondary pathways. These parameters have the general 
form ^y2k^m(0)kNfn(0)"' where k^ = kn, A;_, ^2 corresponds to the nucleation process and 
n is related to the monomer dependence of this process: n = ric — 1, where n^ is the 
critical nucleus size for primary nucleation, n = for fragmentation driven growth and 
n = n2, the secondary nucleus size in cases where monomer-dependent secondary nucleation 
is dominant. The dominance of a single combination of the rate constants implies that 
many of the macroscopic system observables will be correlated since they are dependent on 
the same parameter. A striking examples of this behaviour is provided by the very general 
correlation between the lag-time and the maximal growth rate [HI [T5j, which is manifested 
in the present case in Eqs. (15) and (17) as rmax ~ A and riag ~ A~^. 



Interestingly the rate equations describe sigmoidal curves both in the presence and 
in the absence of secondary nucleation processes. For more complex primary nucleation 
pathways [T6t ITT] the polynomial form for the early time solution is maintained, but higher- 
order exponents are obtained. In the absence of secondary processes, however, the lag-phase 



is less marked since the early time rise is a slower polynomial relationship rather than the 
exponential onset characteristic of secondary pathways [TTj. This observation implies that 
the difference between a high-order polynomial and an exponential may not be apparent in 
experimental data in the presence of noise, and therefore a global analysis of the system 
under different conditions is required in order to obtain robust mechanistic information |14]. 

III. CONCLUSION 

In this paper, we have provided results for the time course of nucleated polymerisation 
for systems that are initially in a mixed state and contain both monomeric and fibrillar 
material. These results generalise the classical Oosawa theory that describes the formation 
of biofilaments to cases where an arbitrary amount of pre-formed seed material is present in 
the system. Furthermore, these results represent a reference to which polymerisation driven 
by secondary pathways can be compared. 
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